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Abstract. We analyzed global properties, radial profiles, and 2D maps of the metal abundances and temperature 
in the cool core cluster of galaxies Hydra A using a deep ~ 120 ks XMM-Newton exposure. The best fit among 
the available spectral models is provided by a Gaussian distribution of the emission measure (gdem). We can 
accurately determine abundances for 7 elements in the cluster core with EPIC (O, Si, S, Ar, Ca, Fe, Ni) and 3 
elements (O, Ne, Fe) with RGS. The gdem model gives lower Fe abundances than a single-temperature model. 
Based on this, we explain why simulations show that the best-fit Fe abundance in clusters with intermediate 
temperatures is overestimated. The abundance profiles for Fe, Si, S, but also O are centrally peaked. Combining 
the Hydra A results with 5 other clusters for which detailed chemical abundance studies are available, we find a 
significant decrease in O with radius, while the increase in the O/Fe ratio with radius is small within 0.1 r2oo, where 
the O abundances can be accurately determined, with d(0/Fe)/d(log]^Qr/r2oo) ~ 0.25 ± 0.09. We compare the 
observed abundance ratios with the mixing of various supernova type la and core-collapse yield models in different 
relative amounts. Producing the estimated O, Si, and S peaks in Hydra A requires either the amount of metals 
ejected by stellar winds to be 3-8 times higher than predicted by available models or the initial enrichment by 
core-collapse supernovae in the protocluster phase not to be as well mixed on large scales as previously thought. 
The temperature map shows cooler gas extending in arm-like structures towards the north and south. These 
structures, and especially the northern one, appear to be richer in metals than the ambient medium and spatially 
correlated with the large-scale radio lobes. With different sets of assumptions, we estimate the mass of cool gas, 
which was probably uplifted by buoyant bubbles of relativistic plasma produced by the AGN, to 1.6 — 6.1 x 10^ M0, 
and the energy associated with this uplift to 3.3 — 12.5 x 10^* ergs. The best estimate of the mass of Fe uplifted 
together with the cool gas is 1.7 X 10^ Mq, 15% of the total mass of Fe in the central 0.5' region. 

Key words. X-rays : galaxies : clusters - Galaxies : clusters : individual : Hydra A - Galaxies ; abundances - 
intergalactic medium - cooling flows 

1. Introduction distribution of several elem ents, in particular iron (e.g. 

iDe Grandi fc Molendil lioOlh and other metals produced 
Clusters of galaxies provide a unique environment for el- ^^^^ supernovae (SN la), which led to the conclusion 
emental abundance measurements and for the study of ^^^^ ^^^^^^^ ^^^^^^ ^^^^ probably due to enrichment 
the chemical enrichment history of the Universe, because SN la in the central dominant galaxies (for a review, 
their large potential wells retain all the metals produced see lWerner et allSi and references therein), 
by the member galaxies. Of particular interest are clus- 
ters of galaxies showing a centrally peaked surface bright- Cooling-core clusters are furthermore in the limelight 
ness distribution and a cool core, whose spectra are of- because of the so-called "cooling-flow problem": the high 
ten richer in emission lines because of the lower central surface brightness in the central peak implies a high den- 
temperatures. In addition, these "cool-core" clusters have sity and a short cooling time, however the observed rate 
been shown to exhibit a central peak in the abundance of cooling of the central gas is often orders of magnitude 

below what is ex pected in the absence of a ny heat sources 



now a Chandra fellow at the Kavli Institute for Particle (for a review, see [Peterson fc Fabian||2006D. As a solution 



Astrophysics and Cosmology, Stanford University, 382 Via it has been proposed that active galactic nuclei (AGN) in 
Pueblo Mall, Stanford, CA 94305-4060, USA the central dominant galaxies can provide enough energy 
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to the central gas to balance the cooling, since signs of en- 
ergetic interaction between the AGN radio plasma and the 
intra-cluster medium (ICM) have b een observed in many 
syste ms (for a recent review, see iMcNamara fc Nulsen 
20071) . Recently, the AGN-ICM interaction has also been 
shown, by theoretical m odels (iRebusco et al. lEooi), hy- 
drodyn amic simulations (IRoediger et al.ll2007l) and obser- 
vations ( Simionescu et al.ll2008bf ) to be a main mechanism 
for transporting the metals produced in the central galaxy 
into the ICM. 

Hydra A was one of the first cooling-core clusters in 
which a displacement of X-ray gas in t he center by radio 



for different co mmonly used abunda nce determinations, 
see the review of IWerner et al.l (|2008( ). 



lobes from the central AGN was found (jMcNamara et al 
2OOOI) and it is still one of the most dramatic exam- 



ples of_AGN interaction ijDavid et al.ll200ll : iNulsen et al 
2OO21 120051 ) . In a deeper Chandra observation of the clus- 



ter, a sharp X-ray surface brightness edge was found 
at radii between 4.3 - 6' (200 - 300 kpc), interpreted 
as a shock wave caused by an AG N outburst with a n 
estimated total energy of 10^^ erg (jNulsen et al.M2005( ). 



which is sufficient to balance radiative cooling for at least 
several 10^ yrs. Further features besides the shock and 
the inner cavities are a ^ 60 kpc long filament running 
from the inner cavities outwards, th at may be reminiscent 



of the "X-ray arms" i n M87 (e.g. iFeigelson et all [19871 



Bohringer et all Il995t iBelsole et all I2OOII : iForman et al.l 



2007h . and surface brightness depressions in the outer re- 
gion inside the shock front, w hich coincide with the outer 
radio lobes ( Wise et al ] |2007l). The relatively low tempera- 
ture (Tx ~ 3.1 - 3.7 keV. lDavid et al.ll200ll ) implies X-ray 
spectra rich in emission lines providing good spectroscopic 
diagnostics. 

The goal of the current work is twofold. Firstly, we 
aim to study the abundances, abundance ratios, and ra- 
dial trends for different chemical elements. This can reveal 
clues about the chemical enrichment history of Hydra A 
and, by comparison with previous results from deep point- 
ings of other nearby bright clusters, can broaden our over- 
all understanding of the origin and distribution of met- 
als in clusters of galaxies. Secondly, we investigate a two- 
dimensional metallicity map to determine what infiuence 
an AGN outburst as violent as the one in Hydra A can 
have on the distribution of chemical elements in the ICM. 

Throughout the paper, we assume Hq = 
70kms-i Mpc-i, Qa = 0.73, and Qm = 0.27. At 
the rcdshift of Hydra A (z = 0.0538), 1" corresponds to 
1.05 kpc. Unless otherwise stated, the elemental abun- 
da nces are given with respect to the proto-Solar values 
of iL^ddiillQoi), the errors are at the la level, and 
upper limits at the 2a level. T he rec ent solar abundance 
determinations by Lodderj ( 2003f l give significantly 
lower abunda nces of oxygen, neo n and iron than those 
measured by Anders fc Grevess3 ( 19891 ). Use of these 
new determinations affects only the units with respect 
to which we present the elemental abundances in our 
paper; the actual measured values can be reconstructed 
by multiplication of the given values with the solar 
normalizations. For a compilation of these normalizations 



2. Observation and data analysis 

2.1. EPIC analysis 

Hydra A was first observed with XMM-Newton on 
December 8th, 2000, for 32.6 kiloseconds (ks). A subse- 
quent 123 ks observation was performed on May 11th, 
2007. Since the second observation is significantly deeper 
and large parts of the first observation were affected by 
soft proton fiares, we will focus in this work primarily on 
the second observation. We extracted a lightcurve for each 
of the three detectors separately and excluded the time pe- 
riods in the observation when the count rate deviated from 
the mean by more than 3a in order t o remove flaring from 
soft protons (IPratt fc Arnaud 2002f) . After this cleaning, 
the net effective exposure is ~62 ks for pn, ^--^81 ks for 
MOSl, and - 85 ks for M0S2. We furthermore excluded 
CCD 5 of M0S2 from our analysis due to its anoma- 
lously hi£h_^ii3^2£_y^:L.§2f^ band during the observation 



(see ISnowden et al. 2008 ). For data reduction we used 



the 7.1.0 version of the XMM-Newton Science Analysis 
System (SAS); the standard analysis meth o ds us ing this 
software are described in e.g. Watson et al. ( 200lf ). 

For the background subtraction, we used a combina- 
tion of b lank-sky maps from wh i ch point sources have been 
excised (jRead fc Ponmanll2003t ICarter fc Readll2007l ) and 
closed-filter observations. This is necessary because the 
instrumental background level of XMM-Newton is vari- 
able and increases with time. Both the blank-sky and the 
closed filter observations were transposed to a position in 
the sky corresponding to the orientation of XMM-Newton 
during our observation. We calculated the count rates in 
the hard energy band (10.-12. keV for MOS, 12.-14. keV 
for pn) outside of the field of view (OoFoV) for our obser- 
vation, for the blank-sky maps and for the closed filter ob- 
servations. These count rates are a good indicator for the 
level of instrumental background in each data set, since 
no photons from real X-ray sources should be recorded 
outside the field of view. For each detector we then added 
to the corresponding blank sky background set a fraction 
of the closed filter data designed to compensate for the 
difference between the OoFoV hard-band count rate in 
the observation and in the blank sky data. We note that 
simply scaling up the blank sky data to match the ob- 
served OoFoV hard-band count rate would implicitly also 
scale up the cosmic X-ray background (CXB) component 
contained in these blank sky maps, leading to an overesti- 
mation of the total backgro und. We compared the blank 
sky background spect ra of Read fc PonmarJ (2003) and 
ICarter fc ReadI (2007) with the current observation in a 
10-14' annulus. We find a good agreement both in the soft 
(0.35-1 keV) and hard (5-10 keV) bands, while between 
1 and 5 keV the spectrum from the observation shows 
an excess consistent with residual cluster emission with a 
temperature of around 2 keV. The best agreement with 
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the 10-14' spect ra, both in th e soft and hard bands, is 
provided by the Read &: Ponman (2 003h backgrou nd for 
the MOS detector and the ICarter fc IleadI (|2007t ) back- 
ground for the pn detector. Consequently, we used these 
respective blank sky fields to estimate our background. 

Out-of-time events were subtracted from the PN data 
using the standard SAS prescription for the extended full 
frame mode. 

2.2. RGS analysis 

We extract the RGS spe c tra fo llowing the method de- 
scribed by IXamura et al. ( 2001 ). We model the back- 
ground using the sta ndard background mod el available in 
SAS (r gsbkgmodel ■ iGonzalez-Riestral 120041) . The cluster 
spectra are extracted from a region which is 3' wide in the 
cross-dispersion direction of the instrument. 

The line emission observed with the RGS from ex- 
tended sources is broadened by the spatial extent of the 
source along the dispersion direction. In order to account 
for the line broadening in the spectral modelling, we con- 
volve the line spread function (Isf ) model with the surface 
brightness profile of the source along the dispersion di- 
rection. We derive the surface brightness profile from the 
EPIC/MOSI image in the 0.8-1.4 keV band. Because the 
radial profile of an ion producing an observed spectral line 
can be different from the radial surface brightness profile 
in the broad band, we multiply the line profile with a scale 
parameter s, which is left as a free parameter in the spec- 
tral fit. The scale parameter s is the ratio of the observed 
Isf width and the width of the Isf model convolved with 
the surface brightness profile. For a flat radial distribution 
of the line emitting ions, the scale parameter is s = 1. 

2.3. Spectral modeling 



Unless otherwise stated, the Galactic absorption col- 



umn density was fixed to A'h — 4.8 x 10 cm 



the 



We use the SPEX package (jKaastra et al.lll996r ) to model 
our spectra with a plasma model in collisional ionization 
equihbrium (MEKAL). We note that the MEKAL model 
implemented in SPEX has a more updated line catalog 
compared to that used by XSPEC; however, SPEX does 
not support any other plasma models, notably APEC. 
Several differential emission measure models which include 
contribution from gas with a range of different tempera- 
tures, rather than a single temperature approximation, 
are available in SPEX. Throughout the paper we will 
mainly use a Gaussian distribution of the emission mea- 
sure (gdem) around the best-fit average temperature, of 
the form 



dx 



^-{x-xq) /2a 



27r(TT 



(1) 



where Yq is the total, integrated emission measure (Yq = 
/ neUYidV), X = logT and xq = logTo with To the average 
temperature of the plasma. Sect. 14.11 gives a detailed ex- 
planation why this is the most appropriate available multi- 
temperature model for the data. 



average value from the two available Hi surveys: the 
Leiden/ Argentine /Bo nn (LAB) Survey of Galactic Hi 



(jKalberla et al.ll2005l = 4.6 8 x 10' 



20 



, ^) and the 

Hi data bv iDickev fc LockmanI (|l990l ) {N-^ = 4.90 x 



10 cm~ ). Point sources identified using the X-ray im- 
ages are excluded from the spectral analysis. The spectra 
obtained by MOSl, M0S2 and pn are fitted simultane- 
ously with their relative normalizations left as free pa- 
rameters. To account for possible gain shifts, the redshift 
was left as a free parameter in the initial fit and fixed to 
its best- fit value when estimating the errors. 

The global spectrum extracted from a circular extrac- 
tion region with a radius of 3' (Sect. 13. ip is fitted in the 
0.35-10 keV band. Because of the low number of counts 
and the high background above 7 keV, the spectra ex- 
tracted from all other regions are fitted in the 0.35-7 keV 
band. The spectra extracted for the analysis of the global 
properties of the cluster and for the radial profiles are 
binned to the optimal bin size using the "obin" command 
of SPEX. This command rebins the data to at least 1/3 of 
the FWHM of the instrument depending on the count rate 
at the given energy. The spectra extracted for construct- 
ing 2D maps of the spectral properties of the cluster have 
much lower statistics and are consequently simply binned 
with a minimum of 30 counts per bin. 

We fitted the normalization, temperature and, where 
the statistics allowed, the abundances of O, Ne, Mg, Si, S, 
Ar, Ca, Fe, and Ni as free parameters. The abundances of 
the other elements heavier than He are fixed to 0.5 of the 
Solar value. 

Due to the high photon statistics of this deep observa- 
tion, our best fit reduced x^s are sensitive to calibration 
problems and to the differences between the individual 
EPIC detectors. To account for this, we include 3% sys- 
tematic errors over the entire energy band used for fitting 
the global properties of the cluster and the radial profiles. 
We usually obtain good (< 2a) agreement between the dif- 
ferent detectors, the only notable exception is mentioned 
in Sect. [3^ 

In addition, to insure an optimal determination of el- 
ement abundances, it is crucial to accurately fit the con- 
tinuum around emission lines. This is not guaranteed by 
available emission measure model distributions, which are 
always a simplification of the true thermal structure in 
any given spectral extraction region. It is necessary there- 
fore, after determining the best average temperature or 
temperature distribution using the full spectral band, to 
fix the thermal structure parameters and fit the elemental 
abundances in narrower bands around the emission lines 
of the respective elements, allowing the spectrum normal- 
ization to vary and thereby correct for small inaccuracies 
in the best determination of the continuum in those nar- 
row energy bands. We accordingly fitted the Si, S, Ar, 
and Ca metallicities in the 1.5-5. keV energy band, and 
the O abundance in the 0.35-1.5 keV band. Unless other- 
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Table 1. Fit results for the EPIC data using a gdem 
model. (Tt is the width of the Gaussian temperature dis- 
tributio n. All abundan ces are relative to the proto-Solar 
units of lLodderd (|2003[ ). 





0-3' 


3 


-8' 


Y (10"" cm"^) 


24.85 ± 0.09 


9.27 


± 


0.05 


kT (keV) 


3.42 ±0.01 


3.65 


± 


0.03 


(7T (logT) 


0.251 ±0.006 


0.210 


± 


0.017 


O/Fe 


0.85 ±0.09 








Si/Fe 


0.65 ± 0.05 


0.45 


± 


0.07 


S/Fe 


0.58 ±0.05 


0.30 


± 


0.10 


Ar/Fe 


0.48 ±0.14 








Ca/Fe 


1.28 ±0.16 


0.47 


± 


0.33 


Fe 


0.445 ± 0.007 


0.324 


± 


0.013 


Ni/Fe 


1.35 ±0.20 








/ d.o.f. 


1449 / 1211 


985 


/ 


972 



wise stated, we report the best-fit Fe abundance using the 
full energy band (both Fe-L and Fe-K lines). 

The Ne abundance could not be constrained with the 
EPIC spectra because the Ne x line lies in the Fe-L com- 
plex where the EPIC cameras cannot resolve the individ- 
ual lines, thus a high Ne abundance can be confused with 
a higher contribution of cooler multi-temperature gas (a 
higher ctt) and vice versa. Consequently, the Ne abun- 
dance was fixed relative to Fe based on the RGS fit (Sect. 
13. 4|) . The Mg abundances could not be determined with 
more than Scr s ignificance in any of the regions. Moreover, 
as reported bv lde Plaa et al. ( 2007 ). the systematic error 
in the effective area calibration of EPIC around the Mg 
energy is so large that the Mg abundance determination 
cannot be trusted. The Ni abundance was not determined 
in any of the regions for which the chosen spectral fitting 
band was 0.35-7 keV, which docs not include the energy of 
the Ni lines. The abundances of Si, S and Fe, on the other 
hand, can be determined with good accuracy in most re- 
gions and, where the surface brightness is sufficiently high, 
we also obtain reliable O, Ar, Ca and Ni abundance de- 
terminations. 

3. Results 

3.1. Global spectra 

To characterise the global properties of the cluster and 
to determine accurate abundance values for as many ele- 
ments as possible, we extract spectra from two large spa- 
tial regions: a circular region with a radius of 3' centered 
on the cluster core and an annulus with an inner radius of 
3' and outer radius of 8'. This should insure good statistics 
in both regions, while keeping the number of source counts 
in the outer annulus well over the background level. 

The best- fit width of the gdem distribution, (Tt, as 
well as all other best- fit parameters are shown in Table [T] 
The ctt's determined from our fits are significantly larger 
than zero both in the inner and outer annuli and the val- 




Energy (keV) 



Fig. 1. The EPIC spectrum of the central 3' and fit resid- 
uals with respect to the best-fit Gaussian emission mea- 
sure distribution model. Data points from pn are shown 
in green, MOS in black, and the best-fit model in red. 



ues are similar to those found bv Ide Plaa et al. 
Sersic 159-03. A spectrum of the central 3' region and 
the best-fit model are shown in Fig. [1] It can be seen 
that the gdem model provides a good fit to both Fe-L 
and Fe-K lines, which a single temperature model did not 
achieve. For a more detailed discussion about the effects 
of the presence of such a multi-temperature structure on 
the spectral properties, see Sect. 14.11 

The deep observation of Hydra A in combination with 
the large collecting area of XMM-Newton make this one of 
the best cluster observations for investigating the chem- 
ical composition of the intra-cluster medium. The spec- 
trum extracted from the circle with the radius of 3' allows 
us to determine relatively accurate abundance values for 7 
elements including Ar and Ca which have relatively small 
equivalent widths (see Table [1]) . We note that there is a 
good agreement between all fit parameters determined us- 
ing the MOS and pn spectra separately: the largest differ- 
ence (2cr) is in determining the Ca abundance. Therefore 
the effect of calibration uncertainties should be small. The 
abundances in the cluster core r elative to Fe are con sis- 
tent with the values measured by[dePiaaeiaO (|2007[) for 
a sample of 22 clusters. The absolute abundance values 
are lower in the outer part of the cluster than in the core. 
The observed relative abundances with respect to Fe are 
also lower in the outer part, although the differences are 
in most cases not significant. 

We do not report the O/Fe value in the outskirts, 
which has a large error and is furthermore uncertain be- 
cause of the low energy wing of the EPIC response which 
does not allow us to resolve the O line well and makes its 
abundance value strongly dependent on calibration uncer- 
tainties and on the uncertainties in the modelling of the 
low energy X-ray foreground, especially in regions of low 
surface brightness. The equivalent widths of the Ar and 
Ca are relatively low and the lower statistics in the outer 
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Fig. 2. Top panel: projected radial temperature profile of 
Hydra A. Bottom panel: radial profiles for four elements 
with the best determined abundance values. The radius at 
which the onset of an abundance p eak is seen is indicated 
by a vertical dotted line. r2oo (from Reiprich fc Bdhringerl 
20021 adapted to Hq = 70kms~i Mpc"!) corresponds to 
23.7'. 



extraction region make the abundance determination for 
these elements uncertain. Except for Si/Fe and S/Fe these 
uncertainties make it difficult to draw strong conclusions 
about the differences in the relative abundances between 
the inner and outer regions. 

3.2. Radial profiles 

We also extracted spectra from 7 circular annuli centered 
on the cluster core in order to determine average radial 
temperature and abundance trends. Our results are sum- 
marized in Table [H 

The upper panel of Fig. [5] shows that the radial tem- 
perature distribution has a dip in the core of the cluster 
and is relatively flat from 1' to 7', with the exception of an 
elevated value in the 3-4' bin, which is lik ely associated 
with the large-scale shock in Hydra A (jNulsen et al..,2002i . 



2 4 
Radius (arcmin) 

Fig. 3. Radial profiles for O obtained by fitting spectra 
from the MOS and pn detectors independently. 



20051 ). For a s tudy of the therma l struct ure associated with 
the shock, see lSimionescu et aL ( 2008a ). Beyond 7', we see 
the onset of a significant temperature decrease. 

The radial distribution of all investigated metal abun- 
dances peaks at the core of the cluster (see lower panel 
of Fig. [2|). Radial gradients are seen in the distribution 
of elements which are predominantly produced both by 
SN la (e.g. iron) and by core-collapse supernovae (SNcc, 
e.g oxygen). Again, the most robustly determined abun- 
dances are those for Fe, Si, and S, but in the inner regions 
(within ~3') where the O abundance determination is less 
sensitive to background subtraction, its value is also quite 
accurate. 

In the last bin, uncertainties in determining the back- 
ground at low energies become important. Thus, we do 
not report the O abundance, which is the most affected 
by such uncertainties, and present in Table[2]the Fe abun- 
dance determined based on Fe-K only. The best-fit Fe 
abundance based on Fe-L and Fc-K combined would have 
been 0.43 ± 0.03 solar (consistent both for MOS and pn), 
which is unlikely at such a large radius from the cluster 
center. 

To confirm the presence of a peak in the distribution 
of the O abundance which, as discussed above, is difficult 
to determine, we also checked the results from the MOS 
and pn detectors separately. The two obtained O abun- 
dance profiles are plotted in Fig. [3l We find a discrepancy 
in the O determination from the two different detectors 
in the most central bin, but beyond this the agreement 
is very good and it is clear that in both detectors the O 
abundances in the inner 3 radial bins (central 2') are sys- 
tematically higher than outside the 2' radius, indicating 
that the peaked O distribution is indeed real. The dis- 
crepancy in the central bin might stem from calibration 
problems at this particular position of the detector, caus- 
ing somewhat different values of ctt (0.19 for MOS and 
0.245 for pn) and different abundances. At other detector 
coordinates, the results agree well. 
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Table 2. Radial profiles obtained by fitting the EPIC data with a gdem model. The emission measure Y was corrected 
for chip gaps and excluded point sources. 





0-0.5' 


0.5-1.0' 


1.0-2.0' 


2.0-3.0' 


3.0-4.0' 


4.0-6.0' 


6.0-8.0' 


Y (10"" cm"'') 


5.15 ± 0.03 


5.82 ± 0.03 


7.25 ± 0.03 


5.56 ± 0.03 


3.66 ± 0.02 


3.74 ± 0.03 


1.89 ± 0.02 


kl (keVJ 


6.17 ± 0.02 


o A A in no 
3.44 ± 0.02 


o A A in no 

3.41 ± 0.02 


3.57 ± 0.03 


3.82 ± 0.04 


3.60 ± 0.04 


3.4d ± 0.08 


(TT (logT) 


0.218 ±0.008 


0.215 ±0.011 


0.246 ± 0.009 


0.258 ±0.013 


0.23 ±0.02 


0.23 ±0.02 


0.08 ± 0.08 


0/Fe 


0.80 ± 0.09 


0.77 ±0.11 


0.85 ±0.13 


0.74 ±0.16 


0.59 ±0.24 


< 0.31 




Si/Fe 


0.60 ± 0.06 


0.67 ±0.06 


0.64 ±0.08 


0.81 ±0.13 


0.52 ±0.18 


0.63 ±0.16 


0.80 ± 0.33 


S/Fe 


0.65 ±0.07 


0.56 ±0.08 


0.51 ±0.10 


0.74 ±0.16 


0.38 ±0.21 


0.53 ±0.22 


< 0.80 


Ar/Fe 


0.85 ±0.18 


0.54 ±0.23 


0.62 ±0.26 


< 0.57 


< 0.69 


< 0.63 


< 0.65 


Ca/Fe 


1.02 ± 0.24 


1.50 ±0.27 


1.82 ±0.31 


1.52 ± 0.49 


1.55 ±0.70 


< 0.68 


< 0.72 


Fe 


0.550 ±0.012 


0.483 ±0.013 


0.392 ± 0.010 


0.308 ±0.012 


0.287 ±0.016 


0.323 ±0.018 


0.29 ±0.07 


/ d.o.f 


1388 / 1222 


1211 / 1107 


1196 / 1067 


1025 / 1029 


972 / 996 


1063 / 1117 


1071 / 917 



David et al.l ( 2001 1) obtain from Chandra data a tem- 
perature profile with a dip in the center down to 3.1 keV, 
increasing up to 3.8-4.0 keV at 3.4', in very good agree- 
ment with what we find in Fig. O Our Fe profile is also 
roughly in agreement with that obtained from Chandra 
data, showing the most pronounced peak in the inner 
~90-150 kpc (1.5-2.5'). However, the XMM-Newton data 
gives much lower Si abundances, and does no t confirm the 
very h igh Si/Fe ratio obtained previously bv iDavid et al 



( 20011 ). While the deprojection used for the Chandra anal- 
ysis could account for relatively higher abundances of both 
Si and Fe compared to our projected radial profiles, it 
would not explain the discrepancy in the Si to Fe abun- 
dance ratios. 



3.3. 2D spectral properties 

To produce two-dimensional maps of the temperature and 
abundance distribution in the cluster, we divided our ob- 
servation into spatial bins with a fixed minimum number 
of counts employing an adaptive binning method base d 
on weighted Voronoi tessellations ( Diehl fc Statler 20061 ). 
which is a generalization of the algorithm presented in 



Cappellari &: CopinI (|2003l ). The advantage of this algo- 



rithm is that it produces smoothly varying bin shapes 
that are geometrically unbiased and do not introduce 
artificially- looking structures. We created a background- 
subtracted count map of the observation in the energy 
range 0.4-7.0 keV, combining all three EPIC detectors, 
and binned this map to 50^ counts per bin for generating 
the temperature map and 150^ counts per bin for gener- 
ating an abundance map, which requires higher statistics. 
The binning to 50^ counts was constrai ned to follow the 
large-scale radio lobes seen at 327 MHz (jLane et al.ll2004 , 
see white contours in the upper panel of Fig. while the 
150^-count bins were constrained to follow the tempera- 
ture contours shown in black in the upper panel of Fig. 
[4l We furthermore constrained the bins in the central re- 
gions (with a high count rate) not to be smaller than the 
extent of the XMM point-spread function (PSF). Since 
the statistics do not allow a reliable determination of any 



abundance except that for Fe, we fixed the O, Ne, Mg, Si, 
S, Ar, Ca and Ni abundances in the fit relative to Fe by 
the ratios representative of the global properties in the in- 
ner 3', as presented in Sect. 13. ll Also, in the bins with 50^ 
counts, the statistics do not allow us to fit the spectra with 
free ctt- We therefore fixed it to 0.2, a typical value found 
in the radial profiles in Sect. 13.21 We note that fitting the 
bins for the temperature map with a single temperature 
model, both with free and fixed iVn, shows the same struc- 
tures. In the bins with 150^ counts used to determine the 
metallicity map, we can obtain best-fit values for ctt, how- 
ever these are very uncertain, resulting in a very noisy ctt 
map, and almost always in agreement with 0.2 (only one 
bin out of 80 is more than 5cr away from 0.2 and three bins 
more than 2.5cr). The obtained metallicity maps with ctt 
free and fixed to 0.2 are very similar and show the same 
features. We therefore choose to present the Fe abundance 
map obtained using ctt = 0.2 in the lower panel of Fig. |4l 
The temperature map clearly shows cooler gas extend- 
ing in arm-like structures towards the north and south. 
The cool gas structure to the north follows well the north- 
ern radio lobe, which rises from the cluster core towards 
the NE and then either bends towards the NW (however, 
the bend in the cool gas feature is stronger than that of the 
radio lobe) or extends out into an umbrella/mushroom- 
shape at the position where the N radio lobe bends. The 
cool gas towards the south is also most probably associ- 
ated with the AGN activity, although the feature in the 
temperature map is oriented at a slightly more easterly 
angle than the southern radio lobe. The Fe abundance 
map also shows features extending towards the north and 
south, in the same direction as the cool gas features, 
suggesting that the cool gas is more abundant than the 
surrounding halo. The cool "arm" to the NE is associ- 
ated w ith a bright 1' long filament found bv lNulsen et al 
(|2005r ) using Chandra data. 



3.4. High resolution spectra 

The spectrum from an extraction region which is 3' wide in 
the cross dispersion direction and effectively ~10' long in 
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Fig. 4. Top panel: temperature map of the Hydra A 
cluster, using a minimum of 50^ counts per spatial bin 
(colorbar units are keV). Temperature contou rs are over- 
plotte d in black, 327 MHz radio contours from lLane et al.l 
are over-plotted in white. The cross marks the sur- 
face brightness peak associated with the cluster center. 
Bottom panel: Fe abundance map, using a minimum of 
150^ counts per spatial bin. Co ntours are over -plotted in 
black (colorbar units are solar, Lodderd 20031 ). tempera- 
ture contours are over-plotted in white. Note that the Fe 
map is a zoom-in of the region shown in the temperature 
map. 



the dispersion direction allows us to determine more accu- 
rate projected 0/Fe and Ne/Fe abundances than possible 
with EPIC. RGS also allows us to put better constraints 



o 
o 




20 

Wavelength (A) 

Fig. 5. Combined RGS I''* and 2"'' order spectrum of 
Hydra A extracted from a 3' wide strip centered on the 
core of the cluster. The continuous line represents the fit- 
ted model. The FeXVll-Fexx lines emitted by cooling gas 
between the Nex and Oviii lines are not visible in the 
spectrum. 

Table 3. RGS fit results using a 3' wide extraction region 
centered on the cluster core. 





2T 


gdem+lT 


Yi {W" cm-^) 


21.58 ±0.22 


21.29 ±0.25 


kTi (keV) 


2.37 ±0.06 


3.42 (fixed) 


(TTi (logT) 


0.0 (fixed) 


0.19 ±0.03 


Y2 (lO"^*^ cm-3) 


0.30 ± 0.06 


0.26 ± 0.04 


kTa (keV) 


0.64 ± 0.04 


0.62 ± 0.04 


0/Fe 


0.74 ±0.10 


0.76 ±0.11 


Ne/Fe 


0.73 ±0.18 


0.84 ± 0.20 


Fe 


0.23 ±0.02 


0.35 ± 0.03 


Scale s 


0.95 ±0.14 


1.01 ±0.20 


/ d.o.f. 


1145 / 915 


1273 / 915 



on the presence of cooling gas in the cluster core than 
EPIC by looking at individual lines of the Fe-L complex. 

The projected global temperature determined by RGS 
is ~1 keV lower than that determined by EPIC. Similar 
differences between the mean temperatures determined 
by EPIC and RG S were found in othe r clusters (e.g. 2A 



0335+096, M 87: 1 Werner et all l2006bl |a[). In the case of 



Hydra A, the reason for this discrepancy could be the fact 
that the RGS fit is performed in the soft-energy band, 
which is not sensitive to emission from higher-temperature 
gas present in the gdem model (see Sect. l4.1.2l for a more 
detailed discussion). To check this, we also performed a 
fit to the RGS data fixing the temperature to the best-fit 
EPIC temperature in the central 3' and leaving the ctt 
free. We obtain a value for ctt comparable to the typical 
results from EPIC. The elemental abundances in this case 
are considerably higher than for the single-temperature 
case, and thus in better agreement with the absolute val- 
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ues determined with EPIC. The abundance ratios, on the 
other hand, are unchanged. The reduced is somewhat 
worse for the gdem fit than for the single temperature fit, 
which is partly due to the low-temperature wing of the 
gdem approximation not being a perfect description of the 
temperature structure and partly due to the fact that the 
RGS/EPIC cross-calibration is not fine-tuned enough and 
the RGS and EPIC spatial extraction regions are differ- 
ent, such that fixing the temperature to the best-fit EPIC 
temperature causes a poorer fit. 

While in several other cooling c ores observed with RG S 
(M 87, 2A 0335+0 96, Centaurus: [Werner et all l2006al lbl: 



Sanders et al the spectral lines from the intermedi- 



ate ionization states of Fe (Fexvii, Fexviii, Fexix, Fexx) 
are well resolved, we do not observe these lines in Hydra 
A. However, the spectrum fitted with a single tempera- 
ture model or with a gdem model leaves strong residuals 
between ~13 A and ~19 A. These residuals can be well 
fitted, both for the single-temperature and gdem models, 
by adding a cool gas component with kT « 0.6 keV and 
emission measure of «1% of the emission measure of the 
hot gas. The presence of this cool gas is detected with 5.0 
and 6.5a significance for the single-temperature and gdem 
models, respectively. The fit results shown in Tableland 
Fig. [5] assume that the abundances of the cooler gas are 
the same as the abundances of the hot gas. 

The best fit O /Fe ratio is consistent with the value de- 
termined by EPIC. At the 3 keV temperature of Hydra A, 
carbon and nitrogen are almost completely ionised and 
their lines cannot be detected by the RGS. 

4. Discussion 

4.1. Spectral effects of the multi-temperature structure 

Throughout this paper, we used a Gaussian distribution of 
the emission measure around the best-fit average temper- 
ature to model our spectra. We will explain in this section 
why this is necessary, what the implications on the best-fit 
spectral parameters are, and we will evaluate how accurate 
the gdem approximation is. 

Usually, the strongest evidence of the presence of 
multi-temperature structure is the shape of the Fe-L com- 
plex, which is located at an energy where the effective area 
of the detectors is high, and whose shape changes drasti- 
cally with the temperature of the emittin g gas, especially 



for t emperatures below ^2 keV (e.g. iBohringer et al 
12004 ^. In Hydra A, however, the shape of the Fe-L complex 
in the EPIC spectrum does not show strong indications 
of the presence of co oler gas as observed in other cool- 
ing core clu sters (e.g . Werner et al. 2006lj : de P laa et al.' 
I2OO&: Simionescu et al... .2008b) . The evidence of a multi- 
temperature structure in this case lies primarily in the 
fact that, when fitting spectra with high statistics, it be- 
comes evident that the best-fit temperatures obtained fit- 
ting independently the soft (0.35-2 keV) and hard (2-7 
kcV) energy bands are in disagreement. Moreover, when 
using the entire energy band, the spectral model cannot 



simultaneously fit both Fe-L and Fe-K lines adequately. 
In particular, the Fe-K line in the data shows a shoulder 
towards higher energies, indicating Fe-K/3 emission which 
suggests the presence of gas at higher temperatures than 
the fitted value from the single temperature model. 

In Fig. [S] we illustrate the difference between the soft 
and hard-band fits for the central 3' spectrum. The up- 
per panel shows the best-fit single-temperature model in 
the soft band, which provides a very good fit to the Fe- 
L complex indicating the lack of a significant amount of 
cooler gas. The lower panel shows this model extrapo- 
lated to the hard band. The clear discrepancy indicates 
the fact that the best-fit single temperature model from 
the soft band does not account for the presence of gas 
at higher temperatures present in the spectrum. For the 
central 3' spectrum, the best-fit soft-band temperature 
is 2.81 ± 0.05 keV, while the hard-band temperature is 
3.86 ± 0.03 keV, approximately the values spanned by 
the temperature profile and the 2D maps (Sect. 13.21 and 
13. 3p . A Gaussian emission measure distribution is one of 
the simplest models which can account for such a mixing 
of emission from plasma spanning a wider range of tem- 
peratures. Using this model significantly improves the fit 
(x^/d.o.f. 1449/1211) with respect to a single temperature 
model (xVd.o.f. 2005/1212). 

4.1.1. Intrinsic multi-temperature structure 

The multi-temperature structure can be partly due to the 
mixing of different regions with slightly different temper- 
atures in the 2D map and partly due to intrinsic multi- 
temperature structure in each of these regions. To disen- 
tangle these two effects, we took the average temperatures 
and the emission measures in all the Voronoi bins with 50^ 
counts within a radius of 3' and, assuming the gas in each 
spatial bin to be single-phase, constructed a total dY/dT 
curve, shown in Fig[7l We over-plot also the best-fit gdem 
model for the central 3' spectrum. As the figure shows, the 
shape of the dY/ dT curve due solely to spatial variations 
in the different bins does also have a Gaussian shape, how- 
ever it is much narrower than the best-fit value from the 
gdem model (full-width at half maximum of about 1 keV 
compared to the 4.2 keV corresponding to ctt = 0.251). 
With the statistical errors on ctt from Table [1] the dif- 
ference is highly significant (~ 35cr). Note that hydrody- 
namic simulations of clusters which have recently under- 
gone energ etic events show a similarly wide dY/dT curve 
(jRasia et al. 2006) as our best-fit gdem model. 

Assuming, instead of isothermality, a gdem distribu- 
tion with (Tt = 0.2 in each spatial region, we find a 
significantly better agreement between the total dY/dT 
curve constructed from the bins in the temperature map 
and the best-fit gdem model for the central 3' spectrum. 
Remaining small disagreements could be explained for ex- 
ample by deviations of ctt from 0.2 in different regions. 
However, Fig [7] does clearly suggest the need for intrinsic 
multi-temperature structure in each spatial bin. 
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According to the best-fit gdem model, roughly 9.5% of 
the total emission measure in the central 3' comes from gas 
between 0.5-1.5 keV, 10.7% from gas between 1.5-2.0 keV, 
and 9.2% from gas between 6-8 keV. In part, this intrinsic 
multi-temperature structure can be due to projection: on 
one hand, the large-scale shock in Hydra A can produce 
the hard-tail in th e dY/dT distribution (for more quanti- 
tative analysis, see ISimionescu et al.ll2008al ): on the other 
hand, cool gas in the cluster outskirts can add to the tail at 
low temperatures. The low density in the outskirts makes 
it however unlikely that all the cool gas arises from pro- 
jection alone and suggests the possible presence of denser, 
unresolved cool gas blobs. To make this statement more 
quantitative: the observed temperature profile in Hydra A 
peaks at roughly 3.8 keV (in agree ment with the depro- 
jected profile of .David et al.l (j200l[ ) from Chandra data). 
An average temperature of 2.0 keV thus is a factor 0.53 
less than the peak temperature. If we co mpare for exam- 
ple wi th the scaled temperature profiles o f IVikhlinin et al.l 
( 2006h . we do not expect the temperature in the outskirts 
to drop to 0.53 of the peak temperature until around rsoo, 
where little cluster emission is expected. Since we consider 
projection along a cyllinder of 3' radius, the volume does 
not increase significantly at large radii {dV converges to 
7r(3')^dr for large r and is larger for smaller r because of 
the curvature of the spherical shell with radius r contained 
inside the considered cyllinder). Thus, unless we have rea- 
son to believe in a very small dT/dr in the outskirts com- 
pared to the cluster center, which scaled temperature pro- 
files do not suggest, dY/dT = ricnudV/dT oc n^dr/dT 
should decrease roughly as the square of the electron den- 
sity beyond ^3^ where projection effects come into play. 
According to lVikhhnin et al.l ( 2006[ ). the density drops by 
a factor of 15 between 0.23r5oo (3') and 7-500, thus the emis- 
sion measure should decrease by at least a factor of ^200. 
Clearly, Fig. [7] shows that the emission measure around 
2 keV is only a factor of a few, as opposed to a factor of a 
few hundred, below the emission measure from tempera- 
tures expected around 3' radius. This shows that projec- 
tion cannot have an important contribution and that the 
observed cool gas in the multi-temperature distribution 
must come from unresolved cool gas blobs in the central 
parts. 



4.1.2. Comparison with other available 
multi-temperature models 

The gdem model provides a significant improvement over a 
single temperature model, yet the real temperature distri- 
bution is probably even more complex, therefore we should 
test the reliability of the Gaussian approximation. 

Unfortunately, other available models except single- 
temperature and gdem are either poorly constrained 
by the current data set or give unphysical results. A 
power-law shaped emission measure distribution (wdem), 
which provides a good fit to other coohng-core clus- 
ter spectra (e.g. Werner et al. 2006b; de Plaa et al.ll2006t 




Energy (keV) 




Energy (keV) 

Fig. 6. Top panel: . EPIC spectrum from the central 3' 
region in the soft band and best-fit single temperature 
model. Bottom panel: EPIC spectrum from the central 
3' region in the hard band and the extrapolation of the 
best-fit single temperature model from the soft band. Data 
points from pn are shown in green, MOS in black, and the 
best-fit model in red. The clear residuals show the need 
for a multi-temperature model. 



Simionescu et al. 2008bf l. requires a very fiat, poorly con- 
strained slope, which is what would be expected if the 
best approximation of the emission measure is Gaussian 
and we try to approximate it with a power-law. 

A two-temperature fit to the central 3' spectrum re- 
quires kTi = 2mto °ot keV and kTa = 5.08^^ ?!^ keV with 



normalizations Yi = 9.8 



+0.8 
-0.5 



and Y2 



-0.13 



cm . This fit gives only a small improvement with re- 
spect to the gdem fit (x^/d.o.f.=1379/1209 compared to 
1449/1211). The best-fit kTi and kTa are approximately 
the minimum and maximum temperatures represented in 
the gdem model within Ictt and the normalizations of the 
two components are similar, which is again what we would 
expect if the best approximation of the emission measure 
is Gaussian and we try to fit it with a two-temperature 
model (in the two-temperature model, the contribution of 
any phase between the two best-fit temperatures can al- 
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Fig. 7. Total dY/dT curve constructed from the bins in 
the central 3' of the temperature map, showing the multi- 
temperature structure: in black, assuming single-phase in 
each bin; in green, assuming intrinsic gdem emission mea- 
sure distribution with (Tt = 0.2 in each bin. The red curve 
shows the best-fit gdem model for the 3' region (Table [T|) . 
In blue and magenta, different best-fit polynomial dem 
models for the 3' spectrum (Sect. |4TT2|) . 



find a very broad multi-temperature distribution, similar 
to that required by the Gaussian model. The 10*'^ degree 
polynomial agrees very well with the gdem model between 
~2-5 keV and falls somewhat faster at lower and higher 
temperatures. The 9**^ degree polynomial agrees with the 
gdem model at large temperatures but shows a hint of 
a bimodal dY/dT distribution with peaks at 2.5 and 6 
keV (however, still with significant amounts of emission 
between 2.5 and 6.0 keV, excluding a simple 2T approx- 
imation). While keeping in mind these small differences, 
we can conclude that the broad best-fit Gaussian emis- 
sion measure distribution is a good simple approximation 
to the best-fit results of the polynomial dem model. 

It is interesting to further note the fact that the dY/ dT 
curves for both the gdem fit to the central 3' spectrum and 
for the 9*"^ and lO"' degree polynomial dem models peak 
around 2.3-2.5 keV, the best-fit temperature from RGS. 
This suggests that the discrepancy between the best-fit 
RGS and EPIC temperatures is partly due to the more 
limited RGS band, which corresponds to low energies and 
is therefore insensitive to the presence of the hard tail in 
the real temperature distribution. 

4.1.3. Confusion with soft excess 



ways be reproduced by combining these two phases with 
appropriate normalization ratios). The gdem model thus 
provides an equally good fit and a more physical interpre- 
tation of the spectrum than the two temperature model, 
because it is difficult to imagine a scenario which would 
generate ~ 2 and ~ 5 keV gas in roughly equal amounts, 
no gas at intermediate temperatures and, implicitly, no 
mixing between the two gas phases. If the 5 keV com- 
ponent comes from the projected shocked gas layer, for 
instance, which would fulfill the requirement for the ab- 
sence of mixing, it should certainly contain much less gas 
than the entire cluster center which would constitute the 
cooler phase. 

The last test is to fit a generic differential emission 
measure model (dem) available in SPEX. For a grid of 
temperatures spaced by 0.1 logT between a user-defined 
minimum and maximum value, this model constructs a 
spectrum corresponding to each grid point, and then can 
apply several methods to determine how to optimally com- 
bine these spect ra in order to reprod uce the data (for fur- 
ther details, see iKaastra et al" I ll996l) . We chose the poly- 
nomial method, which finds the best dY/dT in the form 
of a polynomial of a desired degree. Other methods gave 
either very similar results to the polynomial method or 
were computationally unstable. The disadvantage of the 
dem model is that it cannot fit any other parameters ex- 
cept those related to the dY/dT behavior. We thus had to 
fix the elemental abundances to the best-fit gdem values. 

We show in Fig. [7] the best-fit dY/ dT curves obtained 
from the dem model using 9'^ and 10'^ degree polynomi- 
als. It is immediately clear that, also using this method, we 



Another way to account for temperature differences in the 
soft- and hard-bands (Fig. [B]) is to allow the Galactic ab- 
sorption column density to vary in the spectral fits. A 
best-fit TVh lower than the real value, allowing for emission 
from gas at low energies, in combination with a best-fit 
temperature slightly higher than the real value, can com- 
bine to give a good description of both the soft and hard 
part of the spectrum. If the multi-temperature structure 
is such as we describe above, the single-temperature fit 
with free TVh will be a significant improvement compared 
to the single-temperature fit with fixed TVh- 

We show in Fig. [8] the TVh profile determined fitting a 
single temperature model to the full (0.35-7 keV) energy 
band. The A^h has a dip on the cluster center, which can 
be mistaken for the possible presence of soft excess gener- 
ated either by non-thermal inverse Compton scattering of 
cosmic-microwave photons on relativistic electrons in the 
cluster or by the presence of a soft (~ 0.2 keV) thermal 



component (for a review on soft excess, see iDurret et al 
2008). However, the A^h profile determined fitting a single 
temperature model only to the soft (0.35-2 keV) energy 
band is flat and in good agreement with the value deter- 
mined from H I surveys, disproving this scenario. 

The Gaussian emission measure distribution model 
with Ah — 4.8 x 10^" cm~^ provides a better fit than a sin- 
gle temperature model with free A'h (x^/d.o.f. 1449/1211 
compared to 1718/1211 for the central 3' region). Fitting 
a gdem model with free Ah further improves the fit to 
xVd-o.f.=1291/1210, with a best-fit A^h = 3.92 ± 0.07 x 
lO^*' cm^^, still significantly lower than the Galactic value 
but closer to it than the best-fit A^h using a single- 
temperature fit [imtofi X 10^0 cm-2). Since the best-fit 
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Fig. 8. A'h profile determined fitting a single-temperature 
model to the full (0.35-7 keV, black triangles) and 
soft (0.35-2 keV, red squares) en ergy bands. The val- 
ues determined from H i surveys (jKalberla et al ] 120051 : 
Dickey fc Lockman 1990[ ). are over-plotted with a hori- 
zontal dotted and dash-dotted line, respectively. 



Nn using only the soft band is in agreement with the 
Galactic value (4.74 ±0.07 x 10^" cm^^)^ ^^j^g i^^^^j. 
the full band must reflect some residual fitting problems. 
The most important effect is that the fitting tries to com- 
pensate for the residuals around the 0-edge (around 0.5 
keV, see Fig. [1]) while constraining the average temper- 
ature by the continuum at higher energies; on the other 
hand, there might be some uncertainties in the calibration 
of the spectral slope of the EPIC detectors at low/high en- 
ergies. Alternatively, the low A'h could be due to the pres- 
ence of more cool gas in the real multi-temperature dis- 
tribution than accounted for by the best-fit gdem model, 
although Sect. 14.1."^ suggests that this is not the case. 

4.1.4. The "inverse" Fe bias 

The best-fit Fe abundances from fitting the soft- and hard- 
bands individually with single-temperature models are 
both significantly lower than the best-fit Fe abundance de- 
termined from a single-temperature fit of the full spectral 
band (for the central 3' spectrum, 0.37±0.02 solar for the 
soft band, 0.41 ± 0.01 for the hard band, 0.499 ± 0.007 for 
the full band). The discrepancy with the best-fit Fe abun- 
dance from the Gaussian emission measure distribution 
model (0.445 ± 0.007) is smaller. This can be explained 
in the following way. The strength of the Fe-L complex 
is, for the same Fe abundance, higher for a cooler tem- 
perature. Similarly, the strength of the Fe-K complex is 
higher for a hotter temperature. If the model allows for 
some contribution to the spectrum of cooler and hotter 
gas compared to the average, a lower Fe abundance can 
reproduce the same Fe-L and Fe-K strengths for which a 
larger Fe-abundance would be needed if only gas at the 
average temperature is allowed in the model. This is es- 



pecially important in clusters with intermediate average 
temperatures (2-4 keV), where both Fe-L and Fe-K emis- 
sion lines are seen with relatively similar statistics, thus 
neither of the two lines predominantly drives the fit of the 
Fe abundance. 

This effect is opposite to that described by 



Buote fc FabianI (jl998f ). in which not considering the pres- 
ence of cooler gas in the multi-temperature structure leads 
to an underestimation of the real Fe abundance. That ef- 
fect occurs primarily in cool systems (with an average kT 
below about 2 keV) where the Fe abundance is determined 
largely based on the Fe-L complex, since Fe-K emission at 
those temperatures is very weak. 

In our case, on the other hand, neglecting the contri- 
bution of both cooler and hotter gas in a Gaussian-like 
emission measure distribution leads to an overestimation 
of the Fe metallicity if both Fe-L and Fe-K lines are seen 
with good statistics in the spect ra. We are able th us to 



explain exactly the bias found bv lRasia et al. (2008'), who 



analyzed six simulated galaxy clusters processed through 
an X-ray Map Simulator (X-MAS), which allowed to cre- 
ate mock MOSl and M0S2 observations. They showed 
that the Fe abundances determined from fitting the mock 
X-ray spectra is systematically overestimated for inter- 
mediate (2-3 keV) clusters, while for both hot and cold 
systems, where either only the Fe-L or only the Fe-K emis- 
sion lines dominate, the Fe abundance is recovered with 
good accuracy. 

4.2. Radial distribution of the chemical elements and 
comparison with other clusters 

In recent years, abundance profiles for several elements 
were determined for clusters which have deep observations 
with XMM-Newton and Suzaku. In Fig. H we compare 
the radial distributions of the O and Fe abundances in 
Hydra A with those found in other detailed spectroscopi- 
cal analyses of clusters: 2A 0335-1-096 (obtained by fitting 
the "wdem" differential emission measure distribution , 
no O abundances were presented; IWerner et all l2006bl) 



de Plaa et al 



Sersi c 159-03 (using the "gdem" model 
IIOQJI) , M87 (using the 2T mo del within 2' and the MOS 



data; Matsushita et al. 2003h. Centaurus (usin g the 2T 



model and MOS data; [Matsush ita et al .l l2007ar). Forna x 
(using the 2 T model within 4^ iMatsushi ta et al.|[2b07bl ). 
and A 1060 jSato et al.l 2007). The r2oo values were taken 
from iReiprich fc Bohringe 3 (120021 ) and scaled to Hq = 
70 kms~^Mpc~^. We excluded all abundance values de- 
termined with less than 3a significance. The Fe abun- 
dance increases towards the center in all clusters, only 
Centaurus shows a drop of Fe in the innermost bin (this is 
also the case for Pe rseus, not included in our sample, e.g. 



Sanders et al.ll2004l ). The Fe abundance in the Centaurus 



cluster is about a factor of 2 higher than in the other clus- 
ters, indicating a longer enrichment time scale. 

In Fig. [ini we compare the radial distributions of the 
relative metallicities of O, Si, and S with respect to Fe in 
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Fig. 9. Top panel: The radial distribution of the 
Fe abundance. Bottom panel: The radial distribu- 
tion of the O abundance. The Hydra A data are 
from t his paper, 2A 0335+096 
(20063), Sersic 159-03 
M87 



froi 



from Matsushita et al 



from Werner et al 



.de Plaa et all (|2006[) 
( 2003h . Centaurus from 



Matsushita et al. (|2007a), Fornax from Matsushita et al 



were taken from Reiprich fc Bohringed ( 20o3 ) and scaled 
to Hq = 70 km s~^Mpc~^. All the abundances are given 
with r espect to the proto-SoIar abundances by Lodderj 
(2003). 



all the above named data sets. The radial distribution of 
the relative abundance of Si/Fe is roughly constant within 
0.1r2oo and the ratio is between ~0.6 and ~1.1 Solar. 
There seems to be a drop in Si/Fe in several data sets 
at 0.05r2oo, but the uncertainties of the Si/Fe ratio are 
generally large and they do not allow us to make defini- 
tive conclusions. The radial distribution of the S/Fe ratio 
shows a large scatter between the different clusters, and 
there is no obvious trend with radius. 

Interestingly, the observed trend with radius of the 
0/Fe rati o is much less prono unced than previously de- 



scribed by Tamura et al. ( 20041 ). who found a peaked dis- 



tribution for Fe and a flat radial profile for O, lead- 
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Fig. 10. Top panel: The radial distribution of the Si/Fe 
ratio. Middle panel: The radial distribution of the S/Fe 
ratio. Bottom panel: The radial distribution of the 0/Fe 
ratio. 



ing to a marked increase in 0/Fe with radius (up to 
0/Fe of 10 solar). Also, the trend of O with radius is 
clearly not flat, and there are indications of an increase 
in the O abundance towards the center for most clusters 
plotted in Fig. [9] (exceptions are the last data point in 
A1060 and a central drop in Centaurus, where Fe also 
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shows a lower value in the center). In fact, taking into 
account all data points used in the plots and their re- 
spective error bars, we find that the O abundance signif- 
icantly decreases with increasing radius, with a best-fit 
dO/d(log]^o'^/''20o) = —0.48 ± 0.07, while the increase in 
the 0/Fe ratio with radius is only less than 3cr significant, 
d(0/Fe)/d(logior/r2oo) = 0.25 ± 0.09. The Fe abundance 
does however decrease with radius significantly faster than 
the O abundance {dFe / dilog^^r / r2oo) = -0.72 ± 0.04), 
suggesting that the shght trend in 0/Fe is real. 

The interpretation of the flat O and peaked Fe pro- 
files was, in previous publications, an early enrichment by 
SNcc in the protoclustcr phase, which led to the SNcc 
products being well mixed, and a later continuous enrich- 
ment by SN la in the cD galaxy which continued to enrich 
the ICM on time scales much longer than SNcc , creating a 
peak in SN la products at the cluster center. The flat Si/Fe 
ratio found in many clusters contradicted this picture. Si 
is produced by both SN la and SNcc and according to 
this scenario its radial distribu tion should be sh a llower 
than that of Fe. This prompted iFinoguenov et al. I (|2002h 



to propose that there are two types of SN la with differ- 
ent Si/Fe yields enriching the ICM. The authors proposed 
that, while the SN la with higher Si yields and longer delay 
times contribute mostly within the elliptical galaxy, the 
SNe la with lower delay times and lower Si yields would 
dominate the enrichment in the ICM at larger radii. 

However, Fig. O shows that based on deep cluster ob- 
servations with superior signal to noise ratio, there is ev- 
idence that the radial distribution of O is also centrally 
peaked. The O abundance determination with EPIC type 
CCDs is intrinsically very uncertain and sensitive to many 
systematics - some of them already discussed above - like 
background subtraction, calibration of the oxygen edge, 
A^H, temperature determination, insufficient subtraction 
of the Galactic foreground emission which contains O Viii 
emission. This is especially important in cl uster outskirts 



with low photon statistics. The results of iTamura et al 



(120041 ) are based on relatively short observations with 
XMM-Newton and their significances of the O abundance 
determinations in the outer parts of the individual clus- 
ters are typically only around la. We can conclude that 
the more accurate values shown in Figs. [9] and [10] sug- 
gest different O and 0/Fe radial trends than previously 
assumed. 

We must therefore revise the picture that O /Fe signif- 
icantly increases with radius and that the central abun- 
dance peak in clusters is created almost entirely by SNe la 
and look for a mechanism that could create not only a 
peak in Fe and Si but also a peak in O. Sources of chemical 
enrichment in the ICM are SN la, SNcc, and stellar winds. 
Slight increases in the relative abundanc es of O and Mg 
towards the center have been reported bv lBohringer et al 



( 2005h . who associate the presence of these peaks with 
stellar mass loss in the central galaxy. After a short incur- 
sion into what the observed abundance patterns can or 
cannot convey about supernova yield models, we will test 



this scenario for the case of Hydra A by evaluating the 
contributions from each source of chemical enrichment. 



4.3. ICM abundance patterns and supernova yield 
models 

We plot in Fig. [n] the [Fe/Si]=logio(Fe/Si) against the 
[Fe/0] ratios for all data points available and over-plot the 
expected trends obtained by mixing an increasing amount 
of SN la products into a gas whose initial composition 
is according to SNcc model y ields. We use a varie ty of 
available SN l a models (from llwarnoto et al.l 19991 ) and 



SNcc models ("Tsuiim oto et al.Hl995HNomoto et ahfeood : 



Kobayashi et al...2006) . 

The amount of [Fe/Si] produced by different SN la 
models decreases from WDD3 to W7, WDD2 and WDDl. 
This is shown by the four curves in the left panel of Fig. 
TTJ which all assume the same SNcc initial mass function 
(IM F) weighted avera g e (bet ween 0.07 and 50 Mq) yields 
from lKobavashi et all (|2006f) with an initial metallicity of 
the SNcc progenitor of Z=0.004 (0.2 solar). In the right 
panel of the same figure, we also show the influence on the 
expected [Fe/Si] vs. [Fe/0] curve if we assume different 
initial metallicities of the SNcc progenitor but the same 
SN la model (in this case we chose to show this comparison 
for W7, which is most commonly used in galactic chemi- 
cal evolutior i models). The two dash ed line curves use the 
yields from iKobavashi et al.l (|2006h with Z=0 (bottom) 
and Z=0.02 (t op). Using the IMF we ighte d (10-50 Mp,) 
SNcc yields of lTsuiimoto et al" ( 1995 ) and Nomoto et al.l 
(I2OO6I) (for IMF- weight ed values between 10 and 50 Mq, 



see iWerner et al.l 120081 ) would shift the model curves to- 
ward higher [Fe/Si]. 

We find a large diversity in the collected data, which 
do not strongly favor any particular SN la model. Because 
the 0/Fe radial trends are weak, as discussed above, there 
is no clear tendency for data points from larger radii to 
move along the plotted curves towards lower fractional 
contributions from SN la. Rather, most data points are 
clustered at around 80% of Fe coming from SN la, which 
translates to a contribution of 30-40% of SN la by number. 

Two interesting extreme examples are Hydra A, whose 
very low Si abundance places it at the uppermost limit of 
[Fe/Si] which can be reached by any combination of SN la 
and SNcc models considered, and M87 at the opposite 
end. The high Si abundance in M87 places it at the lower 
limit of Fe/Si that can be reproduced by currently avail- 
able supernova models, and the data points here favor 
the WDDl scenario. Most of the observed values for other 
clusters fall between the WDDl and WDD2 model curves. 
An exciting puzzle is that, as Fig. II II shows, a higher ini- 
tial metallicity of the SNcc progenitor tends to shift the 
expected curve upwards to higher [Fe/Si]. Higher initial Z 



resul t s in a higher [Ne/Fe] a bundance ratio ( Nomoto et al 



2006t IKobavashi et allbood) . The RGS observations how- 
ever show exactly the opposite trend: M87 which has the 
lowest [Fe/Si] has the largest [Ne/Fe] (1.40 ± 0.11 using a 
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2T model, IWerner et ahlbOOGal ) while Hydra A with the 
highest [Fe/Si] has the lowest [Ne/Fe] (0.73 ±0.18 using a 
2T model, this work). 

Since supernovae are the most important sources of 
heavy elements, it should be possible to reconstruct the 
composition of any ICM by an appropriate fraction of 
SN la to SNcc contributions. Understanding if the scat- 
ter of the data points in the [Fe/Si] vs. [Fe/0] plot has a 
physical meaning or is simply due to spectral fitting issues 
remains a challenge for future observational work, while 
trying to reproduce the diversity of abundance patterns we 
see can become a further input for supernova modeling. 

4.4. The origin of the metal abundance peal<s in Hydra 
A 



In Sect. 13.21 we pointed out that the abundances of Fe, 
O, Si and S are larger in the core than in the cluster out- 
skirts. We investigate here how these metal peaks can be 
produced. For this, we must first evaluate the mass of each 
metal needed to create the respective peak. This is given 
by: 



M„ 



X]j 7mct AZmot j-^gasj" 



(2) 
(3) 



where j is an index running over all annuli in the radial 
profile which are part of the peak, AZy^otj — Zy^etj — 
■^mot,out is the average metallicity excess between annulus 
j and the value in the outskirts at which the profile flattens 
out, Mgasj- is the gas mass inside annulus j, and 7niot 
is the mass fraction of the cons idered metal in the solar 
atmosphere, as determined from iLodderd ( 20031 ). 

For each element, we assume as a "base" abundance 
Zout (the abundance value outside the peak, to which the 
radial profiles fiattcn out) the value determined in the 3- 
8' region in Sect. 13.11 For O, the value in the outer region 
is not well constrained, therefore we will use the average 
of the two data points outside 2' where the O abundance 
could be determined in Sect. 13.21 Our estimated "base- 
abundances" are thus 0.32 solar for Fe, 0.15 solar for Si, 
0.10 solar for S, and 0.20 solar for O. If there is a contin- 
uing wBak_£adial_trend_ou^^ 3' for any of the elements 
(see iLeccardi fc Molendil 120081 ) . these 3-8' values are an 
upper limit to the real abundances outside the peak. Given 
all the uncertainty sources involved, we show in this sec- 
tion only calculations based on the best-fit abundances 
without giving formal error bars. The presented results 
depend on the exact abundance profiles but also, among 
others, on the exact supernova yield and stellar wind mod- 
eling and on the IMF, star formation history and average 
stellar abundances in the central galaxy. We present es- 
timates obtained using the best-guess assumptions based 
on available models for these parameters. 

From Fig. [2l it can be seen that the abundance peak 
for all elements is clearly visible inside a radius of 2' (126 
kpc), marked by a vertical dotted line in the plot. We thus 
sum Eq. [2]over the first three bins (0.-0.5', 0.5-1', and 1- 



2') in the radial profile to obtain Mmet for each element. 
Mgas can be determined from the spectrum normalization 
assuming a spherical or spherical shell geometry of the 
emitting region and using Eq. [6l The corresponding metal 
masses in the peak are, in units of 10^ Mq, 2.7 for Fe, 1.9 
for Si, 1.0 for S and 18.5 for O. 

We assume that the enrichment of the metal abun- 
dance peak by the central galaxy took place over the last 
10^" years (approximately between redshift z = 2 and 
the present). This puts lower limits on the rates of metal 
production required to produce the peak, because merg- 
ers between z — 2 and now could have dispersed some 
of the metals out towards the cluster outskirts. However, 
many clusters do requi re enrichment times of ~ 10^" years 



(jBohringer et al ] l2004h . which may imply that mergers are 
not very efficient at disrupting cool cores. 



According to ICiotti et al.l (jl99l[ ). the gas mass con 



tributed by stellar winds as a function of time t can be, 
i n a simplified form which presupposes a Salpeter IMF 
(|Salpeterj 119551 ) and single-age passively evolving stellar 



population, approximated as: 

M^t) » 1.5 X 10-"Lb {t/tH)-^-^ , (4) 

where Lb is the present-day blue band luminosity in units 
of and tn is the Hubble time. The blue-band luminos- 
i ty cor responding to the magnitude presented bv lPetersonl 
( 19861 ) for the central dominant galaxy is 9.2 x lO^^Ls© 
within an aperture of 1.65'(108 kpc diameter at the red- 
shift of Hydra A, which likely includes the light emit- 
ted out to several times the galaxy's effective radius). 
Integrating this equation from t{z — 2) to the present 
(t — tn), we obtain a gas mass contribution from stellar 
winds of 33.6 x 10^ AIq. Assuming the average metallici- 
ties of the stellar population to be solar, we estimate that 
the corresponding masses of Fe, Si, S and O ejected with 
the stellar winds are, respectively, 0.47, 0.28, 0.14, and 
2.2 (in units of 10^ Mq). We caution that the assump- 
tion of solar abundances for the stellar winds is a very 
uncertain one. Firstly, the stellar metallicities in the cen- 
ters of central dominant galaxies can reach super-solar 
abundances. Secondly, elliptical galaxies have old stel- 
lar populations, meaning that their stars formed in the 
presence of less SN la ejecta and the stellar composition 
should be closer to the SNqq abundance ratios. However, 
Kobavashi fc Arimotol (Il999f ) find an average [Mg/Fe] ra- 



tio in elliptical gal axies not very different from solar (1.3 
solar in the units of iLodderdliooa ). Moreover, cDs in cool- 
ing core clusters do show anomalo us blue light emis sion 
and have star- forming regions (see McNamaral Il995l . for 
the case of Hydra A). These young stars (which, being 
young, should have relatively strong stellar winds) must 
have incorporated much more SNIa products than the old 
populations of typical ellipticals. 

If the rest of the Fe in the peak is produced entirely by 
SN la, then assuming the yield from the WDD3 model of 
llwamoto et al.l (|l999l ). which best reproduces the unusu- 
ally high Fe/Si ratio in Hydra A, 2.6 x 10® supernovae are 
required. For the W7 model, in turn, 3.0 x 10® supernovae 
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Fig. 11. The observed [Fe/ S i] against [Fe/O] ratios (logarithmic values). Left: Expected trends for different SN la 
models (Iwamoto et al. I ll999l WDDl-3 in gray increasing upwards and W7 in blue), assuming the SNcc yields from 
Kobavashi et al.l (|2006l) with Z=0.004 for the SNcc progenitor. Right: The effect on the expected trends when using 
different available SNcc models (SN la yields were assu med to f ollow the W7 model). Dashed lines indicate the effect 
of using SNcc models with different Z's from .Kobavas hi et al.l (2006) (0. and 0.02, i ncreasing upwards), dotted lines 
show trends using other pubhshed SNcc yields (above. lTsuiimoto et al. 1995t below. [Nomoto et ah 2006^ . 



would be needed. On a time-scale of lO"'^" years, this im- 
phes an average rate of 2.6-3.0 SN la per century. The cur- 
rent r ate of SN la in elliptical galaxies is ( Cappellaro et al.l 
1999h 

Ria = 0.18 ± 0.06 (lOOyr)-i {10^" Lq)''^ , (5) 

which implies 1.66 SN la per century for the blue lu- 
minosity of the central galaxy in Hydra A, a factor 
of 1.6-1.8 too small compared to the required rate. 
There is however recent evidence that the SN la rate 
in cluster ellip ti cals i s larger than in field ellipticals. 



Mannucci et al 



( 2008h find a SN la rate in cluster 
ellipticals of 0.28tom (lOOyr)-i {10^°Lq)~^, in excel- 
lent agreement with the rate required in Hydra A. 
Further more, the rate of SN la may have been larger in 
the past. iRenzini et al. I (|l993h propose a time-dependence 
of {t/tH)~^ with k spanning 1.1 up to 2. Integrating this 
between t{z = 2) and now, the average SN la rate over 
the past 10^" years could be 2 to 4.1 times greater than 
the present one, and could have produced the observed 
amount of Fe in the peak over the assumed enrichment 
time. 

Using the yields from the WDD3 model, 2.6 x 10^ su- 
pernovae would produce 0.41 x 10^ Mq of Si, 0.25 of S 
and 0.15 of O. Adding the contributions from SN la and 
stellar winds therefore, we still lack 16.2 x 10^ Mq of O, 
1.2 X 10* Mq of Si, and 0.6 x 10* Mq of S compared to 
the estimated total masses of these elements in the central 
peak. Even with the W7 SN la model, whi ch produces the 



highes t amount of O among the models of llwamoto et al 
(|l999l). still 15.9 X 10* Mq of O would be unaccounted 
for solely by adding the contributions of SN la and stellar 
winds. 

One possible solution is to consider that the metal in- 
put by stellar winds is larger than we previously assumed. 



either because of super-solar metallicities in the central 
galaxy or a higher stellar mass loss rate than predicted 
by models (due, for example, to a more complex star for- 
mation history than assumed). Keeping the assumption of 
solar ratios, if the stellar winds would produce around 8 
times more of each metal than assumed above, meaning 
that a correspondingly lower number of SN la would be 
needed to create the Fe peak, we could reproduce well the 
observed masses of O, Si and S in the peak both assuming 
W7 and WDD3 SN la yields. Without the assumption of 
solar ratios but keeping the solar Fe abundance in the stel- 
lar winds and the mass loss rate calculated in Eqn.[4]does 
not present a viable alternative, since the Si/Fe and S/Fe 
in the stellar winds would need to be 5.3-5.4 solar, more 
than is reproduced by any SNcc model. As a combination 
of the two possibilities, if the stellar winds would produce 
only 3 times more Fe than assumed above but would have 
Si/Fe and S/Fe of 2 solar and 0/Fe of 2.8 solar, we could 
also reproduce the observations well. 

An alternative for creating the additional O, Si and S 
is to consider a contribution from SNcc to the metal peak. 
Assuming an unchanged stellar mass loss rate and average 
stellar wind metallicity and using the i nitial mass func- 
tion weighted average SNcc yields from iKobavashi et al 



(2006), approximately 8 x 10* core-collapse supernovae 
would be needed, in addition to the contribution of SN la 
and stellar winds, to reproduce the observed O peak. For 
Si and S, only 6 and 3x10* SNcc are needed, respectively. 
This would mean either that there was a significant con- 
tribution by SNcc in the central galaxy over the last 10^" 
years (note that almost twice as many SNcc than SN la 
are required!), or if the initial enrichment by SNcc in the 
protocluster phase was not as well mixed on large scales 
as previously thought, and some peak in the distribution 
of SNcc products existed prior to further enrichment with 
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predominantly SN la eleme nts. Both cases provide a lter- 
natives to the approach of iFinoguenov et alj (|2002l ) for 
explaining why the radial distribution of Si is not shal- 
lower than that of Fe, as would be expected if only SN la 
contributed to the central abundance peak. 

Additionally, according to recent simulations 



( Domainko et al. 20061 ). ram-pressure stripping of 



cluster galaxies also leads to a stronger enrichment of 
the cluster centers compared to the outskirts, and could 
lead to comparable Fe, Si, S and O peaks. The typical 
spatial scale of this enrichment however (~ 1 Mpc) is 
considerably larger than the extent of the peak in Hydra 
A (- 130 kpc). 

4.5. Gas uplift by the AGN and metal transport into 
the ICM 



As we point out already in Section 13.31 there is a remark- 
able correlation between the 327 MHz large-scale radio 
lobes, cool arm-like extensions in the temperature map, 
and metal-rich filaments in the Fe abundance map. This 
suggests that buoyant radio bubbles produced by the AGN 
uplift cool, metal-rich gas from the central parts of the X- 
ray halo, thereby contributing to the transport and distri- 
bution of heavy elements produced in t he central galaxy 
into t he 
2008bh . 



ICM ( Churazov et al.l 120011 : ISimionescu et al 



The metallicity map shows a clear elongation along 
the N-S direction, with filaments coinciding very well 
with low-temperature features in the central parts of 
the temperature map (Fig. [4]). Similar elongations in the 
metallicity map are seen in simulatio ns of AGN -induccd 
metal-transport in galaxy clusters (R oediger et al. 2007). 
Towards the North, the cool, metal-rich gas overlaps well 
with the rising stem of the northern radio-lobe. Towards 
the South, the low-temperature feature is oriented at a 
slightly more easterly angle than the southern radio lobe, 
and is associated with a weaker increase in metallicity than 
in the N. The southern radio lobe is bent and clearly more 
disturbed than the northern lobe, and the same process 
which leads to the deformation of the radio plasma could 
also account for the displacement between the radio emis- 
sion and the cool filament. 



Simionescu et al.l (|2008br ) used the fact that the gas 



uplifted by the AGN in M87 is multi-phase to infer the 
metallicity, mass, origin and chemical enrichment history 
of this gas. Unfortunately, in Hydra A, the EPIC spectra 
of the cool arm-like features in the temperature map do 
not require additional temperature components. This can 
be explained if the amount of cool gas which has a tem- 
perature around and below «1 keV, where the shape of 
the Fe-L complex becomes clearly different and can indi- 
cate multi-phase structure, is very low. Based on the Fc-L 
shape, one of the very few spectral indications for multi- 
temperature accessible with EPIC, gas above «2 keV is 
difficult to distinguish from gas in hotter phases. Indeed, 
as shown in Section[3i4]where we do detect a small amount 



of cool gas below 1 keV, the spectrum normalization of this 
gas is only 1.2-1.4% of the total emission of the hot halo, 
which is impossible to uniquely identify using the EPIC 
cameras. 

Based on the RGS result, we can try to infer some 
further properties of the cool gas. The RGS extraction re- 
gion we used was 3' wide in the cross dispersion direction 
and effectively ~10' long in the dispersion direction. If we 
assume that the cool gas is concentrated in the cluster cen- 
ter, thus within a 3'-diameter circular region which is the 
2D-projection of a sphere with a radius of 1.5' (94.5 kpc), 
then the mass of the cool gas corresponding to its emission 
measure of 2.6 x 10^^ cm^^ would be 1.8 x lO^Af©. The 
mass is determined using: 

Y = ln.n,dV^n.n,V^l.2nlV 

M = '^miUiV « (tohtih + mHc"Hc)V" « l.AmnnnV (6) 

i 

where rio, nn and nuo are the electron, proton and helium 
number densities, respectively, mjj is the proton mass, 
muc = 4rnH, and V is the volume of the emitting region. 
Approximately, ric ~ 1.2 riu, and assuming a solar He/H 
fraction, nHc/"-H ~ 9.5%. We neglect the mass contribu- 
tion of ions heavier than He which are scarce compared to 
H and He. 

If we assume on the other hand that the cool gas comes 
only from the cool "arms" seen in the maps and not from 
the entire 1.5'-radius sphere, we can approximate the vol- 
ume of this gas by a sum of cylinders for the Northern 
and Southern "arms" , depicted in Fig. [121 Both cylinders 
are bound to lay inside the 3 by 10' RGS extraction re- 
gion. The N cylinder has a radius of 0.35' and a length of 
1.3', the S cylinder has a radius of 0.3' and is 2.5' long. 
This yields a volume of 8.9 x 10^^ cm'^, and a mass of 
5.1 X 10^° Mq. This is only a rough estimate because it 
assumes the cold gas to be distributed uniformly in the 
considered cylindrical regions. 

We furthermore need to take into consideration the 
possibility that the cool gas may not fill the entire assumed 
volume of the "arms" , but may be concentrated in thin fil- 
aments which occupy only a fraction (so-called filling fac- 
tor) of this volume. Since these filaments are unresolved, 
the only way to estimate this filling factor is to assume 
pressure equilibrium between the cool and hot phases, in 
the form UckTc = rihfcTh which, using n^h oc \/i^,h/14,hj 
implies 



(7) 



where c and h are subscripts denoting the properties of the 
cool and hot gas determined by RGS, respectively. From 
Eqn. m we obtain V^/Vh = 1.0 x lO^^ ^ Fc/Kot, which 
leads to a mass of cool gas of 1.6 x 10^ A/©. This is a factor 
of ^3 larg er than the mass of uplift ed cool gas calculated 
in M87 bv lSimionescu et all (|2008bl ). 

If we assume that the average displacement of the mass 
of cool gas determined above is on the order of 1', which 
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Fig. 12. Zoom- in on the temperature map in Fig 2] 
smoothed with a 16" Gaussian filter. A circle with radius 
of 1.5' and the 3 by 10' RGS extraction region are overlaid 
in white. The cluster center is marked with a black cross 
and the black rectangles mark, in projection, the cylinders 
by which we approximate the volume of the cool "arms" . 
Colorbar units are keV. 



is about 63 kpc, and use the integrated mass profile pre- 
sented bv lPavid et al. ( 2001 ). we can estimate the gravi- 



tational energy needed to upHft this gas at 3.3 x 10^° ergs, 
which is 8 times greater than the energy needed for gas up- 
lift in M87 but nevertheless very small compared to other 
processes related to AGN-ICM interaction in Hydra A, 
which require approximately 10^^ ergs (larg e-scale shock, 
Nulsen et al."2005|; IS imionescu et al.ll2008aL and cavities. 
Wise et al. 2007). 

Finally, we estimate the amount of Fe which is being 
transported by the AGN through the uplift of the cool 
gas. Assuming a gas mass for the cool gas of 1.6 x 10^ Mq 
and a metallicity, as determined from RGS, of 0.35 solar, 
and using M-pc = 7Fc-^^gas, this gives Mpo = 7.9 x lO^M©. 
However, this is almost certainly too low, since the up- 
lifted cool gas should be more metal-rich than the ambi- 
ent gas, as the structure in the Fe map suggests, while 
the abundances of the hot and cool gas were coupled in 
the RGS fit. Moreover, some of the uplifted metal-rich gas 
may have higher temperatures, therefore its mass is not 
included in our estimate of 1.6 x 10^ Mq of uplifted gas, 
which is based solely on the RGS best-fit normalization of 
the 0.62 keV gas. This leads to an underestimation of the 
mass of uplifted gas and implicitly of the mass of uplifted 
Fe. 

Another way to estimate the mass of transported Fe 
is to assume a typical metallicity outside the metal-rich 
regions and apply Eq. [2] summing over all bins which lie 
inside the rectangles by which we approximated the pro- 
jected shape of the N and S arms in Fig. [121 Assuming 
the average metallicity outside the metal-rich regions to 
be 0.445 solar (the average in the inner 3', Sect. 13. 1| ). the 
mass of uphfted Fe would be 1.7 x 10^ Mq, while if the 



Table 4. Uplifted gas mass, Fe mass, and energy, *' with 
equal metallicities of the 0.62 keV and the hot gas from the 
RGS fit, **) if the average metallicity of the total uplifted 
gas is 2 solar. 





gas mass {Mq) 


Fe mass [Mq] 


energy (ergs) 


0.62 keV 


1.6 X 10** 


7.9 X 10^*^ 


3.3 X 10^** 


total 


6.1 X 10^**' 


1.7 X 10^ 


1.25 X 10^^ 



metallicity outside the metal-rich regions is 0.32 solar (the 
"base" abundance for Fe discussed in Sect. 14. 4| ). we obtain 
an upper limit of 6.9 x 10^ Mq for the mass of Fe trans- 
ported by the AGN together with the uplifted gas. 

In turn, the total mass of Fe in the central 0.5' re- 
gion is 11.5 X 10^ M0, meaning that the AGN transported 
about 15% and up to 60% of the amount of Fe currently 
present in the inner part of the Hydra A halo out to larger 
radii by uplifting central, metal-rich gas. We note that the 
0/Fe value determined by extracting a spectrum in the N 
arm (0.80 ± 0.14) is consistent with the average 0/Fe ra- 
tio i n the inner 3' (T a ble (T]). This means, in agreement 
with lSimionescu et al.l ( 2008bh . that the AGN, at least at 
this advanced stage of the evolution of the galaxy, trans- 
ports metals into the ICM without altering the relative 
abundance patterns. 

We already noted that some of the uplifted gas may 
have temperatures higher than 0.62 keV, therefore its mass 
is not included in our estimate of 1.6 x 10^ Mq, which re- 
sults in underestimating the uplifted gas mass. We showed 
a method of estimating the uplifted mass of Fe indepen- 
dent of the estimated mass of the 0.62 keV component. 
If we then assume an average metallicit y for the uplifted 



gas of 2 solar (as in the case of M87, iSimionescu et al 



2008bl ). we can compute the total mass of uplifted gas 
to be 6.1 X 10^ Mq, corresponding to an uplift energy of 
1.25 X 10^^ ergs. The results for different methods of es- 
timating the gas and Fe mass and the uplift energy are 
summarized in Table [J] 

5. Conclusions 

We analyzed a deep ~ 120 ks XMM-Newton exposure of 
the cooling core cluster of galaxies Hydra A. We extracted 
spectra from two large regions in the cluster core and in 
the outskirts to study the global properties of the cluster 
with the best statistics possible. We also analyzed the RGS 
spectrum extracted from a 3 by 10' region. We find that 

— the shape of the Fe-L complex in the EPIC spectrum 
does not show strong indications for the presence of 
cooler gas as observed in other cooling core clusters. 
However, a multi-temperature model is needed to si- 
multaneously fit both Fe-L and Fe-K lines appropri- 
ately. The best available model achieving this is a 
Gaussian distribution of the emission measure (gdem) 
around the best-fit average temperature. 
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— the best-fit {gdem) model shows a very broad distribu- 
tion, with a full-width at half-maximum of 4.2 keV. We 
also fitted polynomial emission measure distributions 
which confirm the broad shape of dY/dT. The dis- 
tribution is significantly broader than expected from 
the substructures visible in the 2D map in the consid- 
ered region, suggesting the presence of intrinsic multi- 
temperature structure in each bin of the temperature 
map. Wc suggest that this multiphase structure can be 
due in part to the projection of shocked gas in front 
of and behind the cluster center, in part to the projec- 
tion of cooler gas from the cluster outskirts along the 
line of sight, and in part to dense, unresolved cool gas 
blobs in the considered extraction region. 

— wc can accurately determine abundances for 7 elements 
in the cluster core with EPIC (0,Si,S,Ar,Ca,Fe,Ni). In 
the cluster outskirts, only Si, S and Fe abundances are 
determined with better than 3cr significance and the 
large errors do not enable us to draw conclusions about 
possible differences in Si/Fe and S/Fe with respect to 
the cluster center. 

— the abundances of 3 elements (0,Ne,Fe), one of which 
is inaccessible with EPIC, are determined with RGS. 
The 0/Fc ratios from EPIC and RGS are consis- 
tent. While no multi-temperature structure beyond the 
gdem model can be constrained with EPIC, the RGS 
fit requires a cool component with a temperature of 
0.62 ± 0.04 keV and a normalization of only 1.2% of 
the hot ambient. This gas is detected with a signifi- 
cance of 6.5(7. 

— the Gaussian emission measure distribution model 
gives lower Fe abundances than a single temperature 
model, which can explain why simulations show that 
the best-fit Fe abundance in clusters with intermediate 
temperatures is over-estimated. 

We also determined temperature and abundance profiles 
from seven annuli and compared our results with radial 
profiles of other clusters for which deep observations and 
detailed chemical enrichment studies are available. We 
show that 

— the abundance profiles for Fe, Si, S, but also O are 
centrally peaked in Hydra A. 

— the radial profiles of the O abundance are peaked also 
in other clusters for which deep data are available. The 
increase in O /Fe with radius is very small. Combining 
the Hydra A results with 5 other clusters for which 
a detailed chemical abundance study has been per- 
formed, we find dO/d(logior/r2oo) = -0.48 ± 0.07, 
while the increase in the O /Fe ratio with radius is only 
less than 3a significant, d{0/Fe)/d{logiQr/r2oo) = 
0.25 ±0.09. 

— stellar winds and the chemical enrichment by the num- 
ber of type la supernovae needed to produce the Fe 
peak in Hydra A do not reproduce the estimated O, 
Si and S peaks. For this, either the amount of metals 
produced by stellar winds would have to be 3-8 times 
higher than predicted by available models or 3 — 8 x 10^ 



SNcc would be needed in addition to the contribution 
from stellar mass loss and SN la. Possibly, the initial 
enrichment by SNcc in the protocluster phase was not 
as well mixed on large scales as previously thought, and 
some peak in the distribution of SNcc products existed 
prior to further enrichment with predominantly SN la 
elements. 

— mainly because of the low Si abundance. Hydra A re- 
quires either a WDD3 or W7 SN la model to reproduce 
the observed relative abundance patterns. Most other 
clusters lie between the WDDl and WDD2 models. A 
30-40% contribution by SN la compared to SNcc (by 
number) is needed to reproduce the observed relative 
abundance patterns at all radii (below 0.1 r2oo) in all 
the considered clusters. 

— the best- fit Galactic absorption column density, A^h, 
when fitting the full energy band, is lower than the 
value determined from H l data and has a minimum on 
the cluster center, which may be interpreted as inverse 
Compton or soft excess emission. The best-fit A'h in 
the 0.35-2 keV band however is in agreement with the 
H I data, dismissing this possibility. 

Finally, we produced 2D temperature and metallicity 
maps of Hydra A. From these we can conclude that 

— the temperature map shows cooler gas extending in 
arm-like structures towards the north and south. The 
cool gas structures, and especially the northern one, 
appear to be richer in metals than the ambient medium 
and spatially correlated with the large-scale radio 
lobes. The northern "arm" is associated with a bright 
1' long filament seen in the Chandra image. 

— based on the geometry of the cool "arms" seen in the 
temperature map and on the best-fit normalization of 
the 0.62 keV component in RGS, the estimated mass 
of cool gas, which was probably uplifted by the AGN 
to create these structures, is 1.6 x 10^ Mq. The en- 
ergy needed for this uplift is 3.3 x 10^* ergs, which is 
8 times greater than the energy needed for gas uplift 
in M87 but nevertheless very small compared to other 
processes related to AGN-ICM interaction in Hydra A 
(large-scale shock and cavities), which require approx- 
imately 10^^ ergs. 

— the best estimate of the mass of Fe uplifted together 
with the cool gas is 1.7 x 10^ Mq, 15% of the total 
mass of Fe in the central 0.5' region. If the average 
metallicity of the uplifted gas is 2 solar, as in M87, 
the total mass of uplifted gas (not only the 0.62 keV 
component) is 6.1 x 10^ M© and the uplift energy is 
1.25 X 10^^ ergs. 

— the 0/Fe value in the N arm is consistent with the 
average 0/Fe ratio in the inner 3'. The transport of 
metals by the AGN thus presently does not alter the 
relative abundance patterns. 
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